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^ ' Abstract 

o ; 

^^ ■ We analyse a classical model of gravitation coupled to a self interacting scalar field. 

f^ , We show that, within the context of this model for Robertson- Walker cosmologies, there 
exist solutions in the spatially non-flat cases exhibiting transitions from a Euclidean to 

^^ ■ a Lorentzian spacetime. We then discuss the conditions under which these signature 

■*««« , changing solutions to Einstein's field equations exist. In particular, we find that an 

qh| upper bound for the cosmological constant exists and that close to the signature changing 

J ' hypersurface, both the scale factor and the scalar field have to be constant. Moreover we 

bi). find that the signature changing solutions do not exist when the scalar field is massless. 

> 

%. 

1 Introduction 

The principle of causality is one of the most fundamental pillars on which classical as well 
as quantum physics is based. It exhibits itself in terms of the causal structure of a four- 
dimensional Lorentzian space-time manifold. Theories like quantum fields rely on this structure 
for a sensible interpretation of their results in a fixed background geometry. However, there 
is no apriori reason which would suggest that this causal structure should remain unaffected 
by the dynamical equations of general relativity. If we accept that in quantum cosmology [|l| , 
amplitudes for gravity can be expressed as sum of histories of 3-geometries @] with different 
global 4-topological structures, then this may give rise to structures that cannot have pure 
Lorentzian or Euclidean geometries [Q and this implies a possible signature transition [Q. For 
example, an interesting feature of present day quantum cosmology is characterised by its 
use of Riemannian signature spaces to arrive at an observable Lorentzian-signature spacetime. 
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Therefore, in order to gain further isight, it might be useful to study the dynamics of this 
phenomenon in simple models. 

Traditionally, one of the features of classical gravity is that the signature of the metric is 
usually considered as fixed. This is not a property demanded by the field equations, but rather 
a condition one imposes on the metric before looking for solutions to Einstein's equations. If 
one relaxes this condition, one may find solutions to the field equations, which when interpreted 
suitably, exhibit a signature transition. 

In general, for studying signature dynamics in the classical realm of gravity 0, one assumes 
that there exists a hypersurface which separates the space into two disjoint regions, one being 
Euclidean and the other Lorentzian. On this hypersurface, henceforth called the junction, the 
metric is either continuous and necessarily degenerate, or discontinuous with nondegenerate 
one-sided limits [^. Therefore one can categorize this problem according to the continuity 
of the metric. Having chosen the form of the metric, there are two approaches for solving 
the field equations. The first (continuous) approach is to assume the validity of Einstein's 
field equations throughout and identify the form of the field equations at the junction with 
the junction conditions 0, H, 3. In the second (discontinuous) approach one finds solutions 
by solving Einstein's equations in disjoint regions next to the junction and then uses specific 
junction conditions {i.e. Darmois junction conditions) to match the solutions |jlO|- This 
furnishes an alternative way of categorizing the problem. 

These different approaches have resulted in mainly two different sets of junction condi- 



tions which have led to a controversy between various authors [^. Within the framework 
of Robertson- Walker cosmology using either approach, the scale factor, the scalar field and 
their first time derivatives are required to be continuous across the junction. However within 
the continuous approach, Hayward using a continuous metric has obtained extra conditions 
restricting the first time derivatives to be zero at the junction. The controversy mentioned 
above stems from the claims by some authors that these extra junction conditions are unnec- 



essary, see e.g. |T2|, [T^ and therefore the matter conservation need not hold under signature 
change |]13| . 



Recently, Kossowski and Kriele [14], within the continuous approach, have proved a lo- 
cal existence and uniqueness theorem for solutions of Einstein's equations with dust energy 
momentum tensor in a class of m-dimensional signature changing spacetimes where initial con- 
ditions are given on the hypersurface of signature change. They have also proved a similar 
theorem in the case where the energy-momentum tensor represents a scalar field which is not 
self-interacting. 

In this paper, we study a classical model where signature transition is possible, using the 
continuous approach with a continuous metric. This model is the same as that proposed 
by Dereli and Tucker ^ in which a real massive scalar field is taken as the matter source 
interacting with gravity and itself in a Robertson- Walker geometry whose signature evolution is 
controlled by a preferred coordinate. One then seeks solutions which are smooth and continuous 
across the junction on which the metric is degenerate. For the spatially flat universes, this 
results in exactly solvable Einstein's field equations |^. Here we show that signature changing 
solutions also exist for arbitrary curvature in Robertson- Walker cosmologies. We do this by 
finding analytic solutions to the field equations close to the junction, and confirm and extend 
them away from the junction by a numerical method. We also show systematically the general 



set of restrictions that need to be imposed on the junction conditions and the parameters of 
this model including the cosmological constant. The general set of allowed conditions and 
parameters thus obtained include the ones used in reference [0 as a special case. In particular 
we find an upper bound on the cosmological constant, and the need for the scalar field to be 



massive in order to have continuous solutions. The present work is also complementary to ||T^ 
in that it solves the problem for a self-interacting scalar field. Having chosen a continuous 
approach, we obtain a set of junction conditions which, not surprisingly, are more restrictive 
than those obtained in the discontinuous approach. However, the restrictiveness of the junction 



conditions obtained here lies somewhere in between those obtained in P| and [|I4[. We shall 
return to this point in the conclusions. 

2 Field equations 

We begin by briefly reviewing the model, and for comparison purposes, use the same notation 
as that of 0. Consider gravity coupling to a scalar field through 

G = kT[0], (1) 

where the scalar field is a solution of 

dU 
A^-- = 0. (2) 

Here, G = Ric — ^glZ is the Einstein tensor constructed from torsion-free connections com- 
patible with the metric, and U is a scalar potential for the real scalar field interacting with 
itself and gravity through the stress-energy tensor T[0], given by 

T[0] = d(P0d(P- ^g{d(P, d(P)g - U{(P)g. (3) 

The above coupled equations are to be solved in a domain that would lead to Robertson- Walker 
cosmologies with Lorentzian signature. However, if the metric is suitably parametrized, one 
can expect to see continuous transition to a Euclidean domain. As in |^, we adopt a chart 
with coordinate functions {P , x^ , x'^ , x^} where the hypersurface of signature change would be 
characterised by /3 = 0. The metric can be parametrized to take the form 

, = -Pdp 0dp+ ^^ f^^fly,^, E dx^ ® dx\ (4) 

where r^ = X^^^iX*, and depending on the values oi k = { — 1,0, 1} one would have an open, 
flat or closed universe, respectively. Now, it is apparent that the sign of /? determines the 
geometry, being Lorentzian if /? > and Euclidean if /? < 0. For /? > 0, the traditional cosmic 
time can be recovered by the substitution t = (2/3)/?^/^. Adopting the chart {t,x*} and using 
equations (|lD through (^ with units in which k = 1, one finds 

= TT + W), (5) 




fR\ fkV k (^2 

'[r)^'[r) ^W - -| + W). (6) 

where a dot represents differentiation with respect to t. These equations are not all inde- 
pendent. Equation (H) can be obtained by combining equations (|]) and (|^). Equations (^ 
through d^) can also be derived in the context of Lagrangian dynamics. One defines the action 
4-form 



A 



^n-^g{d<p,d<p)-ui<p) 



* 1, (8) 



whose variations with respect to g and (p give equations (|Tp and (|^), respectively. This action 
can be written as 

A = CdtAdx^ Adx^ Adx^, (9) 

where 

Cdt = [-3RR^ + 3kR + i?3[0V2 - f/(0)]} rft + d{3R^R). (10) 

After dropping the total derivative term and using the following transformation, 

X = R^/'^cosh{a(j)), (11) 

Y = R^/^smh{a(l)), (12) 

we obtain a mechanical analogue for this system given by the Lagrangian 

2a''Cdt = j-X^ + y2 + ^(x2 _ y2)i/3 _ 2^,2(^2 _ Y^)U{(I){X, Y))\ dt, 



(13) 



with — oo<(/)<cx3,0<-R<oo and a^ = |. By a straightforward Legendre transformation 
we obtain the "Hamiltonian" 



2a^ndt = j-X^ + r^ _ ^(x2 _ y2)i/3 ^ 2a\X'' - Y^)Ui<P{X, Y))\ dt, 



(14) 



which we identify with the total energy density of the system. If we transform this expression 
back to the form of R and (p, we see that equation (|^) is nothing more than a "zero energy 
condition." Of course any solution of the Euler-Lagrange equations would yield a constant 
total energy. However, Einstein's equations demand zero energy solutions only. 
The above Lagrangian is considerably simplified if we take the potential to be 

2a^{X^ - r2)f/(0(x, Y)) = aiX^ + agF^ + 2bXY, (15) 

where ai, 02 and b are free parameters. This potential is required to have natural characteristics 
for small 0, so that we may identify the coefficient of 0^/2 in its Taylor expansion as a positive 
m^ and U{(f) = 0) as a cosmological constant A. In terms of the above parameters, A = U\(j,=o = 
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ai/(2a^) and m^ = d'^U j d(t?\^=Q = ai + 02- The features of this potential have been discussed 
in more detail in 0. 

The dynamical equations for X and Y in terms of the evolution variable (3 now become 



Y" = \ (^^ r - \pkY{X' - Y^r^'^ - [5{a,Y + 6X), 
X" = ^ (^ X' - \(3kX{X^ - r2)-2/3 + p^^^x + bY), 



(16) 
(17) 



subject to the "zero energy condition" 



2a^n = ( ^ J (-X'' + r ') - -k{X^ - F2)1/3 + ^aiX^ + a2Y^ + 2bXY) = 0. (18) 

Here, a prime represents differentiation with respect to (3. The coupled equations ([T6|) and 
(|17|) must now be solved and, since their solutions render H a constant in /?, equation (|18D 
becomes only a restriction on the initial conditions. In the next section we first find analytic 
solutions to these equations which are valid close to the junction, and then find the solutions 
for the full range of /5 by a numerical method. 

3 Solutions 

3.1 Analytic solutions close to the junction 

In this section we find analytic solutions which are valid near the junction, and obtain the 



restrictions imposed by equations (p^)-([T8|) on the junction conditions. This is done by noting 
that in order to have well behaved solutions close to (3 = 0, the first term of equation (piSf ) 
shows that we must either have X'{j3) ~ (3^'' and Y'{i3) ~ /?"'', where n^, Uy > 1/2, or 
|X'(0)| = |y(0)|. However, the first terms on the right hand side of equations ( |T6D and (|l^ 
impose a more severe restriction. These two equations admit solutions X'{P) ~ /?^/^ and 
Y'{P) ~ /5^/^ close to /5 = 0, however, these class of solutions do not admit real or C^ solutions 
across /3 = 0. One can show that regular solutions close to /? = are of the form 



X{(3) = A,l3^ + Xo where A, = l 

y 



° +a^Xo + bYo 



Y{/3) = Ay(3' + Yo where A 



y 



4 (X^ - ro2)2/3 

3 kYp 

'^{Xl-Y,y/^ 



02^0 - bXi 



(19) 
(20) 



with Xq = X{0) etc. Therefore, the initial conditions on the first and second derivatives must 
satisfy the relations 

X'(0) = F'(0) = and X"{0) = F"(0) = 0. (21) 

Strictly speaking the conditions on the second derivatives are not initial conditions but 
rather consistency checks, since we have coupled second order equations. The above relations 



confine the possible solutions to a restricted class. Therefore, the initial values for the functions 
X and Y must now satisfy, c.f. equation (O) 



- h{Xl - Y^Y'^ + {a^Xl + a^Y^ + 2hX^Y^) 



0. 



(22) 



Although equation (|2^) is equivalent to a sixth order algebraic equation which cannot be 
directly solved analytically, we can solve it by going back to the original variables R and 0. 
The solutions are either -R(O) = giving 0(0) = ±oo, which we exclude because we have been 
seeking continuous solutions across /5 = 0, or -R(O) ^ (it is a free parameter) with 



0(0) 



2a 



cosh 



DB ± bVW^^WT¥ 

52-62 



(23) 



where 



D 



9k 



fll — 02 



AR{oy 



and B 



ai + a2 



m 

T 



The contour plots of equation ( P^ for A; = ±1 are given in figure 1. Along the contours, one 
finds the possible initial values for X and Y . Of course, the acceptable values of X and Y ( 
X > \Y\) can also be obtained analytically from equation (PB]). 

In order for this model to have a valid quantum extension, at least in the scalar sector, we 
need to require the potential to have a minimum so that a real mass can be defined. That is 



d^U 



m„ 



min 


= m^i 


A- 


462 


where 


(Pmin ~ 


1 , 1 26 
:---tanh ^. 
2a m^ 



Therefore we need to require rn^ > 2|6|. This requirement along with the requirement of having 
real values of X{0) and Y{0) satisfying the zero energy condition and leading to real values 
for i?(0)(> 0) and 0(0) give the following restriction on the parameters 



9A; 



ai 



AR{Oy 



<\b\< 



fli + a2 



It is important to note here that equations ( p!9ip0| ) imply 

X(0) = and r(0) = 
X(0) : 



^A^ and YiO) = ^Ay. 



These in turn lead to 

m = 
m = 



and 0(0) = 0, 

— — and 0(0) = -^sinh[2a0(O)] cosh[2a0(O)]. 

R{0) 2a a 



(24) 

(25) 
(26) 

(27) 
(28) 



The physical significance of these constraints and their relevance to the present controversy on 
the junction conditions will be discussed in the last section. 



3.2 Numerical solution 

A noticeable feature of equations (|T6|) and ( ]T7| ) is that ioi k ^ they are singular at the 
critical values of f3 (/3c) at which Y{/3c) = ±X(/3c). In terms of the original variables we have 
R{Pc) = and 0(/?c) = ±C)cQ. Any attempt in solving these equations involves handling these 
"moving singularities," as one encounters them when integrating the coupled equations. To 
handle these singularities, we establish jump conditions across them as follows: close to /5 = Pc, 
we assume that the solutions have the following linear forms 

X± = a± + b±(3, 

Y± = c± + d±/3, 

where ± refers to the right or left hand sides of the singularity respectively. Substituting the 
above equations in ([T6| ) and (|l^) and integrating in the interval (3c — ^, Pc + £> where 2e is the 
distance across the jump, gives the following equation at Yc = ±Xc 

9 (2X fV^^ 
h-b- = -l \l^^_y,A = ^(d^ - ^-)' (29) 

where e can be taken as small a value as is desired for any required accuracy. Moreover, the 
values of a± and c± are determined by the requirement of the continuity of X and Y at (3c- 
Equation (|29| ) establishes our jump condition for handling the singularities of the differential 
equations. For ease of comparison with the solutions of the spatially flat {k = 0) Robertson- 
Walker universe, we use the same set of parameters as in ||^ (6 = 2, A = 0, m^ = 4.5) and 



choose our initial conditions according to equations (0) and (|^). We note that since equation 
(|I^) is a constant of motion, if it is satisfied aX (3 = 0, it will be satisfied at all other values of/?. 
For integrating equations ([T6|) and (p^Tf), we have used the fourth order Runge-Kutta method. 
The resulting solutions for k = ±1 are shown in figure 2. As a measure of the accuracy of 
the solutions, we have computed the "zero energy condition," equation (]T8|), as a function of 
(3 for k = ±1. We have found that the values of "total energy" stay very close to zero in the 
full range of (3 shown in figure 2, thus indicating the validity of our numerical solutions. As a 
further check, we have numerically recovered the analytic solutions presented in for k = 0. 
In figure 3 the variations of (p and R and in figure 4 that of the curvature scalar are shown as 
a function of /5. 

4 Conclusions 

We have shown that in the context of classical general relativity, using a model proposed in 
0, there could exist signature transitions in Robertson- Walker cosmologies with arbitrary 
curvature. In this model, these transitions are embodied in the form of smooth continuous 
functions describing solutions to Einstein field equations coupled to a scalar field. These solu- 
tions describe a geometry in which the covariant metric tensor is degenerate on a hypersurface 
characterised hj (3 = and undergoes a transition from a Euclidean to a Lorentzian signature. 



^None of the values of Pc coincide with the junction at which /3 = 0. 



As in the k = case, the singular behaviour of the scalar field characterizes the beginning of 
the Euclidean domain. This feature is also refiected in the behaviour of the curvature scalar. 
For k = the solutions for X and Y are not singular anywhere in the Lorentzian domain. 
However, for k ^ 0, there are singularities a.t Y = ±X but they are very mild and show up 
only at the second derivative level. As k increases we observe more rapid oscillations in X and 
Y, leading to more rapid variations in all the physical quantities. However, for all values of k, 
they vary extremely slowly in the vicinity of the signature changing hypersurface, as can be 
seen from figures 3 and 4. 

As far as the controversy on the junction conditions in the current literature is concerned, 
we find that in our continuous approach to signature change (with a continuous metric) the 
restrictions on the initial conditions are given by equations ( pT| , |27|p8D . These conditions are 
slightly more restrictive than those of [Q , because equation ( PT| ) implies that the second deriva- 
tives of R and with respect to the evolution parameter (/3) are zero, whereas the analysis 
presented there does not seem to require this. However, the junction conditions of |T^ seem to 



be too restrictive since it requires the second (cosmic) time derivatives of R and to be zero, 
contrary to equation (^). As was mentioned before, we only require the second derivatives 
with respect to (3 to be zero. 

We now discuss the physical significance of the restrictions on the parameters of the model 
which is presented in equation (p^. When this equation is written in terms of the physical 
parameters it becomes 

3A 9k ,,, rri^ 

T - 4iP(0) - '"' - -■ ""' 

In particular we have Aixiax = 3A;/-R^(0) + 2?7i^/3 and this implies that A cannot have an 
arbitrarily large positive value. Rather more importantly, note that as m ^ 0, equation (|30D 
implies that 6 = and subsequently equation ( p3D shows that 0(0) —>■ ±oo. Therefore, in order 
to have continuous solutions across the junction, the scalar field has to be massive. We expect 
this result to be independent of the particular choice of the potentual U{(f)). 

It is interesting to note that although the solutions for k = and k = ±1 are qualitatively 
not very different, there is a pronounced difference in the curvature scalar. That is 7^(0) = 
Qk/R\0). 
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Figure Captions 



Figure 1: The contour plots of the allowed values of X and Y, satisfying the equation of 
constraint (^2]) at /3 = for A; = ±1 — figures 1(a), 1(b). The point (0,0) is a solution and 
the curves approaching this point actually pass through it, although this is not shown on the 
plots due to the limitations on the numerical accuracy. 

Figure 2: Solutions for X (broken curve) and Y (solid curve) as a function of (3, for k = ±1 
— figures 2(a), 2(b). The values of the parameters are 6 = 2, A = 0, m^ = 4.5. 

Figure 3: Variation of (p and R with (3 for 6 = 2, A = 0, m^ = 4.5, and k = 1 (broken curves) 
and A; = — 1 (solid curves). Note that (p and R vary extremely slowly in the vicinity oi (3 = 0. 

Figure 4: Variation of the scalar curvature 7^ with /3 for 6 = 2, A = 0, m'^ = 4.5, and k = 1 
(broken curve) and k = —1 (solid curve). Note that 7^(0) = 6k/[R'^{0)]. 
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